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Abstract: Two concepts of symmetry for the distributions of positive ran- 
dom variables Y are log-symmetry (symmetry of the distribution of log Y) 
and R-symmetry [7]. In this paper, we characterise the distributions that 
have both properties, which we call doubly symmetric. It turns out that 
doubly symmetric distributions constitute a subset of those distributions 
that are moment-equivalent to the lognormal distribution. They include the 
lognormal, some members of the Berg/Askey class of distributions, and a 
number of others for which we give an explicit construction (based on work 
of A.J. Pakes) and note some properties; Stieltjes classes, however, are not 
doubly symmetric. 
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1. Introduction 

Let Y be a random variable following an absolutely continuous distribution F 
with density / on the positive half- line R + . One natural concept of symmetry of 
this distribution is that Y/6 and 6/Y have the same distribution for some S > 0. 
Since this is equivalent to ordinary symmetry, about log 5, of the distribution of 
logy, we will refer to this as log-symmetry ([6] Section 12. A. c, [ L]). In density 
terms, 

y 2 f{&y) = f (J^j , for some S > and all y > 0. (1) 

Note that S is the median of the log-symmetric distribution F. 

An alternative interesting concept of symmetry on R + is the R-symmetry 
introduced and investigated by Mudholkar and Wang [7] . This is defined directly 
in terms of the density and another constant 8 > by 

f{6y) = f (^j , for some 9 > and all y > 0. (2) 
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In fact, 9 is the mode of F (if / is unimodal). A comparison of these two concepts 
of reciprocal symmetry was given by Jones [5] . 

The focus of this note is on distributions on K + which have both properties 
(1) and (2), log-symmetry and R-symmetry; we will call such distributions dou- 
bly symmetric. It is easy to see that no continuous distribution on K + can be 
doubly symmetric with the same constant (8 = 9). However, as noted in [5], the 
lognormal distribution — with density 

1 {\ogy~n) 2 \ 

2 V 2 J (3) 

in its usual parametrisation resulting from taking the exponential of a random 
variable following a normal distribution with mean fi and variance a 2 — is both 
log-symmetric about 5\ = e M and R-symmetric about 9\ = e^~ a . In that paper, 
the first author of this note was "tempt [ed] to conjecture that the lognormal is 
unique in this respect" — this proves to be wrong! 

In this note, we characterise the class of absolutely continuous distributions 
on R + which are doubly symmetric (call this class DS). It turns out to be a 
proper subset of the wide class of absolutely continuous distributions on R + 
which arc moment-equivalent to the lognormal distribution (e.g. [9], Section 
3; use MEL to denote this class); moreover, this subset of MEL distributions 
contains more densities than just the lognormal. See Section 2 for our main 
result. Our work will rely heavily on deep and detailed work of Pakcs [8, 9]. 

The first, and most famous, MEL subfamily is a version of the Stieltjes class 
with density 

fetv) = fo(y) [1 + "in {2ir(\ogy - fj) /a 2 }] , (4) 

— 1 < e < 1. This is the class considered by Stieltjes [12] and Heyde [4]; see 
[13] for a more general formulation. It is intriguing to find that such Stieltjes 
densities are not DS. In fact, members of the Stieltjes class (4) are neither log- 
symmetric nor R-symmetric except when e = 0. Instead, y 2 f t (5\y) = f- e (5\/y) 
and f c {9iy) = f~ e (9i/y), the former comprising Theorem 3.1 of [!)]. And it is 
clear that no other values of 6 or 9 can be found to make / e log-symmetric or 
R-symmetric, respectively, for e ^= 0. 

For examples of MEL densities that are DS, and some investigation of their 
properties, see Section 3. 

The paper finishes in Section 4 with some further comments. 




2. Main result 

Write k = 5/9. Note that k > 1 because the median is not less than the mode 
for either log-symmetric or R-symmctric densities [5]. Then, manipulating (1) 
and (2) shows that DS densities satisfy 

1 f (v\ = y 2 f(y) (5) 

k 2j \k 2 J 9 2 k 4 { ' 
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(together with cither of (1) and (2)). This, therefore, says that DS densities are 
the same when you rescale by k 2 (left-hand side of (5)) and when you perform 
weighting with weight function y 2 (right-hand side of (5); note that E{Y 2 ) is 
indeed equal to 9 2 k 4 ). 

Rescaling, let g(x) = 0k 2 f{9k 2 x) be the density of X = Y/(6k 2 ). Then, (5) 
translates to 

^■9 (p) = y 2 9(y) (6) 

which is the same as (5) except that E(X 2 ) = 1. Now, (6) is identical to (3.15) 
of [8] provided we set Pakes's q = k~ 2 . Next, write h(z) = g(y / z)/(2^/z) to be 
the density of Z = X 2 . Then, as Pakes notes, the problem of solving (6) for g 
is equivalent to solving 

F*(FH fc M (7) 

with E{Z) = 1. 

Pakes ([8], Theorem 3.1) gives the general form for any distribution function 
solving (7) with E(Z) = 1. These distributions, it turns out [10], are a subset 
of the distributions that are moment-equivalent to the log-normal law. We are 
interested only in absolutely continuous versions of these distributions, for which 
Pakes [8] writes down the form of the density, h* say, at his (3.10) (see also 
Theorem 3.3 of [!)]). It has a piecewise form on intervals (k~ i ^ l+1 \ k~ Al ], i = 
0,±1,±2, .... Explicitly, 

OO 

h*(y) oc ^ k 2t( - l+ ^y l uj{k M y)I{k- i ^ +1 '> <y< fc~ 41 ) (8) 

i— — oo 

where u> is a nonnegative function on (fc~ 4 , 1] such that h* is integrable. Since 
f(y) = {2y/e 2 k A )h (y 2 /{6 2 k A )) , solutions of (5) take the form 

OO 

/*(y)oc £ e" 2l k 2l ^y 2l+1 uj{e~ 2 k^ l ~^y 2 )I{9k~ 21 <y <9k 2 - 2% ). (9) 

i— — oo 

However, densities /* constitute a superset of DS densities since the latter 
also need to satisfy one of (1) or (2). 



Theorem 1. Absolutely continuous doubly symmetric densities have the form 
f* given by ( 9) with u chosen to satisfy 

^>( u ) = ^(^r)> k^<u<i. (10) 



Here, 



1p{u) = UUj(u). 
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Proof. First, write (9) as 

oo 

/*(!/) oc ^ 2i fc 2i(i+i y i+i -^~ 2 fc 4(i ~V) 

i— — oc 

where JF(n) = uj(u)I(k~ A < u < 1). Now, consider requirement (5), recalling 
that <5 = fe0. We have 

OC 

*/ 2 /*(fc%)oc £ fc^V^^" 2 */ 2 ). (11) 

i— — oo 

Also, 

OO 

/*(fc(9/y) oc Yl kW)y-^n^- 2 /y 2 ) 

j=-oo 
oo 

= kW-Qv"- 1 ^*;-"- 2 /y 2 ). (12) 

i— — oo 

Expression (12) follows from the previous line by setting j = —i, the choice 
determined to be the one that gives the T functions in (11) and (12) the same 
support. It follows that (11) = (12) if 

fcV^-¥) ^(fc- 4i - 2 A, 2 ) 

which, by setting u = k 4l ~ 2 y 2 , is equivalent to (10). The reader can verify that 
the same result is obtained by considering (G) in place of (5). □ 

Remarks. Note that the mode 9 enters (9) simply as a scale parameter. Note 
also that, although /* is the density of the (scaled) square root of a random 
variable whose distribution is MEL, it remains the case that /* itself is MEL; 
this is because the (scaled) square root of a lognormal random variable is itself 
lognormal. It will also prove useful, in Section 3, to write 

oo 

f*(9y) oc ^- 1 )y 2 -V(fc 4(l - 1) y 2 )/(fc- 21 < V < k 2 ~ 2i ). (13) 



3. Examples of DS distributions 
3.1. Lognormal distribution 

Lognormal density fo given by (3) or any alternative scaling of it by factor c 
say, has E(Y 2 ) = c 2 e 2 ^+ 2ff2 = 9 2 k A since 9 = ce^-* 2 and k = . Note that 
the i = 1 term of (13) shows that, on support (fc~ 4 ,l], ip(u) oc f*(9y/u)/y/u. 
Then, the log-normal's psi-function is 

^oc^expf-^lM (14) 
y ' ^ I 8 logfc J v 1 

which satisfies (10). 
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3.2. Another known MEL distribution 



For a first example of a class of MEL distributions some of whose members 
are DS and are not the lognormal, we turn to the Askey/Berg class ([1, 2, 8] 
Section 3, especially pp. 834-835, [9] Section 3). One way of writing the density 
of this class, employing a particular scaling for convenience, is as 



l 



/ 7 (y) cx yi-i/L-kiy) where L k (y) = ]T y n k~2 n , (15) 

n— — oo 

7 £ M, y > 0, is Ramanujan's theta function (e.g. [3]). Now, 



L k {k c y) 


OO 


1 2 , 

-ooV K 2 


Lk(k c /y) 


Z_-m— - 


. eo y-"A;-i« a + cw 




Z— /n— - 











1 2 , 



^ 9 1 2 ' 

Here, J^ 1 sums over j £ {. . . , —2 — c, — 1 — c, — c, 1 — c, 2 — c, . . .} and J^ 2 sums 
over j € {■ ■ ■ , —2 + c, — 1 + c, c, 1 + c, 2 + c, . . .}. The two sums are equal if 
c is integer or half-integer, in which case (1) holds with 5 — k 1 and (2) with 
8 = fc 7_1 if 7 is also integer or half-integer. This chimes with (1.10) of Theorem 
1 of [3] whose integer n is our 2c. We have thus shown that Askey/Berg densities 
with integer or half-integer values of 7 are DS. 

In fact, limited computational comparisons of the lognormal and Askey/ Berg 
densities with the same median and mode (as well as moments) have failed to 
find any visible differences between them. (Often, the two densities are identical 
to the eighth decimal place, sometime to just four.) The version of the lognormal 
density that is compared with / 7 at (15) has [i = jlogk, a 2 = logfc and hence 
density 

My) = v^m y exp l"2ioifc 

It follows that, to a very good degree of approximation, it must be that 

1 log 2 y 



L k (y) ccexp ^- bgfc 

Indeed, it can be proven that / 7 and fo are exactly equal at all points of the 
form k p , p = 0, ±1, ±2, .... To see this, note that, exactly, 

J- 2 , 1 2 » , 1 / \2 1 2 

~ p) cx eV . 



Lk(k p ) = J2 fc ~ 2 ™ +np = elP ~ J2 fc ~ 2(i; 



n— — oo 
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3. 3. Constructing further specific DS MEL distributions 

Appropriate functions ip satisfying (10) can be constructed by specifying non- 
negative ip arbitrarily on [fc~ 2 , 1] and extending it uniquely to the remainder of 
[fc~ 4 , 1] by denning 

V>W = V(^), k- 4 <u<k- 2 . (16) 

Note that ip is continuous at k~ 2 . Also, it can readily be shown from (9) that 
continuity of density /* at the "grid endpoints" k 2 \ i = 0, ±1 ± 2, . . . requires 
0(1) = ?/>(fc~ 4 ) — equivalent to Lemma 3.1 of [9] — which is immediate from 

(10) ; it follows that if ip is chosen to be continuous then /* is continuous. We note 
Pakes's ([9], p. 1278) observation that continuous MEL densities "probably are 
the exception" and that this construction seems a promising way of extending 
this class (with double symmetry thrown in). 

Now let ip be continuously differentiable. From (16), continuous differentiabil- 
ity of /* at "grid midpoints" fc 2 *" 1 , i = 0, ±1, ±2, . . . requires ip'(k~ 2 ) = 0. From 
(13), continuous differentiability of /* at grid endpoints k 2l ,i = 0, ±1,±2, ... 
requires 2ip'(l) = — ip(l) < 0. 

Unimodality of /* with unique mode 8 demands that: 

(i) u l ~^ 1 ' 2 'ip(u) be increasing for all i = 1,2,... which follows if ^/uip>(u) is 
strictly increasing; 

(11) u l ~^/ 2 ^(u) be decreasing for all i = 0, —1, —2, . . . which follows if ili{u)/^/u 
is strictly decreasing. 

It can then be shown that, when ip is differentiable, these two requirements 
reduce to choosing ip on (fc~ 2 , 1) such that |(log-0)'(M)| < (2w) _1 . 

One parametric family of unimodal DS (MEL) distributions arises if we take 

xj)(u) = tp a {u) = u a ~ {1/2 \ k~ 2 <u<l, 0<a<l. (17) 

We note at once that, while continuous, the resulting density is not continuously 
differentiable at both grid mid- and end-points for any value of a. However, this 
at least makes for densities that more obviously differ from the lognormal. To 
this end, we briefly investigate the simplest case, that of a = in (17). The 
corresponding density is 

r( , Y.Z-ook^Ky/e? 1 - 1 ^- 21 <y/e <k 2 - 2 ^) 

Ju(y) e{2\ogk + Y,f =l r 1 k- 2 ^{k^ ' [ ' 

It is pictured in Figures 1 and 2 in conjunction with the moment-equivalent 
lognormal density in each case. Density (18) acts as a piecewise (odd-power) 
polynomial approximation to the lognormal with discontinuities in derivative 
at grid endpoints. Discrepancies between the two densities are clearest at their 
common modes. 
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Fig 2. Densities f^(y) (in black) and fo(y) (in red) for k = 1.75 and 8 = 0.1 (solid lines), 
8 = 0.5 (dashed lines), 8 = 1.0 (dotted lines) and 8 = 2.0 (dot-dashed lines). 
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4. Further remarks 

That doubly symmetric distributions are moment-equivalent to the lognormal 
distribution is essentially immediate from the derived relationship (5) from 
which we get 

r2(s+2) 

E(Y S+2 ) = -0j^jE(Y s ), s£R. (19) 

(Inter alia, the existence of all positive and negative moments of / is thereby 
established.) This relationship is satisfied by the lognormal distribution with 
5 = cxp(/i), 9 = exp(^i — a 2 ) and E(Y S ) = exp(s/x + \s 2 a 2 ). In the closely 
related context of (7), Pakes [8, 9] notes that the ratio of the moment functions 
of / and of the lognormal is periodic. 

Here is another take on periodic differences between solutions. Define h(w) = 
log f(e w ), w £ R. Then, a little manipulation based on (5) yields 

h(w — 2 log A) = 2w- 21og6»- 21ogfc + h(w). (20) 

Thus, h comprises a periodic function, with period 21ogfc, plus a linear one. 
Moreover, for DS densities we also have, from (2), that 

h(w) = h(2log9- w), (21) 

displaying symmetry of h about 2\og9. Combining (21) and (20) yields 

£(w-2\ogk) =£(2log6-w) (22) 

where £(w) = h(w) — in. 

The periodicity evident in (20)-(22) can be utilized to formulate a character- 
ization of the lognormal distribution within the class DS. If X has a lognormal 
distribution then for any 7 > 0, X 1 is also a lognormal variable and conse- 
quently is also DS. If one assumes only that for any 7 > 0, X 1 £ DS, it follows 
that the functions in (20)-(22) have multiple periodicities and, as a consequence, 
the function h(w) must be quadratic. In this manner we may formulate: 

Theorem 2. If X 1 £ DS V7 > 0, then X has a lognormal distribution. 

We conjecture, but are as yet unable to verify, that it will suffice to replace in 
the hypothesis the statement "V7 > 0" by the statement "for two judiciously 
chosen values of 7" . 
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